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ABSTRACT 

The growth rate of long wavelength kinetic instabilities arising due to the interaction of a collimated 
beam of relativistic particles and a cold unmagnetized plasma are calculated in the ultra relativistic 
limit. For sufficiently culminated beams, all long wave-length modes are shown to be Weibel-unstable, 
and a simple analytic expression for their growth rate is derived. For large transverse velocity spreads, 
these modes become stable. An analytic condition for stability is given. These analytic results, which 
generalize earlier ones given in the literature, are shown to be in agreement with numerical solutions 
of the dispersion equation and with the results of novel PIC simulations in which the electro-magnetic 
fields are restricted to a given k-mode. The results may describe the interaction of energetic cosmic 
rays, propagating into the far upstream of a relativistic collisionless shock, with a cold unmagnetized 
upstream. The long wavelength modes considered may be efficient in deflecting particles and could 
be important for diffusive shock acceleration. It is shown that while these modes grow in relativistic 
shocks propagating into electron-positron pair plasmas, they are damped in relativistic shocks propa- 
gating into electron-proton plasmas with moderate Lorenz factors Fgh ^ (frip/'nT'e)^ ■ If these modes 
dominate the deflection of energetic cosmic rays in electron-positron shocks, it is argued that particle 
acceleration is suppressed at shock frame energies that are larger than the downstream thermal energy 
by a factor of > Fgh- 
Subject headings: shock waves acceleration of particles cosmic rays 



1. INTRODUCTION 

Current understanding of gamma-ray burst (GRB) 
"afterglows," the delayed low energy emission follow- 
ing the prompt 7-ray emission, suggests that the ra- 
diation observed is the synchrotron emission of ener- 
getic non-thermal electrons in the downstream of an 
ultra-relativistic collisionless shock driven into the sur- 
rounding interstellar medium (ISM ) or stellar wind 
(|Zhang fc Meszarosl[200l |Piran|[200l . 

This model requires a strong magnetic field and a 
large population of energetic electrons to be present 
in the downstream. Observations suggest that the 
fraction of post-shock thermal energy density car- 
ried by non-thermal electrons, e^ , is large, e^ « 

0.1 (e.g. [Zhang fc Meszaro d I200I iFrail et al.l [200TI : 



role in generating and maintaining the inferred magnetic 
fields. 

The main challenge associated with the downstream 
magnetic field is that the field amplitude must remain 
close to equipartit ion deep into the downstream, over dis- 
tances ^ lO^^lsd (jGruzinov fc Waxmanlll999t iGruzinovl 
l2001a( ). While near equipartition fields on skin depth 
scale are likely to be produced in the vicinity of the shock 
by electromagnetic ( e.g. Weibel-like) instabilities (e.g . 
Blandford fc Eic hled fl98 7l: IGruzinov fc WaxmanI [19991: 



Medvedev fc Lo eb 1999^: iWiersma fc Achterbergri2004[ l . 



iFreedman fc Waxman 200i nBerger et al.l 120031 ). The 
fraction of post-shock thermal energy carried by the mag- 
netic field, eb, is less well constrained by observations. 
However, in cases where es can be reliably constrained 
by multi waveband spectra, values close to equipart ition, 
63 -^ 0.01 to 0.1, are inferred (e.g. lFrail et"al1[2000h . 

The non-thermal energetic electron (and proton) 
population is believed to be produced by the dif- 
fusive (Fermi) shock accel e ration (DSA) mechanism 
(for reviews see iDruryl 119831: iBlan dfordfc Eichle^ 119871: 
IMalkov fc O'C Drurvll200l y 

The required magn etic fields in the shock frame in 
the downstream (e.g . iFrail et all I2000D and upstream 
(|Li fc WaxmanI 1200 6( ) regions are much larger than the 
ambient field, and thus require substantial amplification. 
The accelerated particles are likely to have an important 

|itay.rabinak@weizniann.ac.il| 



they are exp ected to decay within a few skin-depths 
downstream ([Gruzinovl l2001a| ) . This suggests that the 
correlation length of the magnetic field far downstream 
and possibly upstream must be much larger than the skin 
depth, L :^ Isd, p erhaps even of the order of the distance 
from t he shock ( Gruzinov fc WaxmanI 119991 : IGruzinovl 
[200Tat IKatz et al.il2nn7n . 

The search for a self-consistent theory of collision- 
less shocks has led to extensive numerical studies 
using the particl e i n cell (PIC) based algo r ithms (e.g. 
Gruzinovl I2001allbl: IMedvedev et al.l [200a ISilva et ah! 



20031: iNishikawa et al.l 120031 iFrederiksen et all 12004 : 
Jaroschek et all l2004t iSpitkovsk^ [200l [2008a : 
Martins et al.l 120091 ) . Such simulations have pro- 



vided compelling evidence for acceleration of particles 
and generation of long lasting near-equipartition mag- 
netic fields. However numerically simulating the long 
term behavior is challenging and is currently restricted 
to pair (e~^e~) pl asmas in 2D (e.g. iSpitkovskvl l2008bl : 
'Keshe t et aLllBol ). 

Large scale magnetic fields may possibly be gen- 
erated in the upstream by the interaction of the 



beam of CRs propagati ng ahead o f the s hock and the 
upstre am plasma (e.g. iKatz et all 120071 : iKeshet et al.l 
l2009f ). In particular high energy CRs naturally in- 
troduce large scales due to their large Larmor ra- 
dius, and the large distances to which they prop- 
agate into the upstream. Instabilities arising from 
the interaction of rclati vistic beams an d cold plasmas 
have long been studied (jAkhiezeii 119751 and references 
therein) and are suspected of am plifying the magnetic 
field in the shock transition layer (Gruzinov & Waxmad 
1999; Mcdyedcv & Locb 1 999; Wiersma & Ach t crberg j 
200l iBret et all 120051: iLvubarskv fc Eichled 1 2006 : 



Achterberg et al.ri2007l: lAchterberg fc Wiersmal 12007 : 
BretJl2009D . In (jLemoine fc PelletieHl2009[ ) a systematic 



study of these instabilities for the lowest energy CRs, 
with energies comparable to the thermal energy of the 
shocked plasma, and their application for Fermi acceler- 
ation is given. 

In this paper we analyze long wavelength plasma insta- 
bilities resulting from the counter-streaming flow of high 
energy CRs ,7 ^ Fgh, running far ahead of the shock and 
a non-magnetized upstream plasma. The analysis is re- 
stricted to long wavelength modes, k <C wq/c, which are 
expected to deflect particles efficiently. For simplicity it 
is assumed that the particle distribution is homogenous. 
The paper is organized as follows. In § [5] we calculate 
the growth rate of long wavelength modes. We sepa- 
rately discuss highly coUimated beams and beams with 
a significant transverse velocity spread, and derive a con- 
dition for the stability of these modes. In § [3] we discuss 
the possible implications of these results to collisionless 
shocks. In §5] we summarize the main results and conclu- 
sions. An estimate of the saturation level of the modes is 
beyond the scope of this paper. Throughout this paper, 
units with c = 1 are assumed (c is retained in some of 
the expressions). 

2. ANALYSIS 

Consider a homogenous, anisotropic distribution of 
particles consisting of a cold plasma and an axi- 
symmetric beam of ultra relativistic CRs. The analysis 
is carried out in the rest frame of the cold plasma which 
initially has zero magnetic and electric fields. 

The plasma frequencies, in this frame, of the cold 
plasma and the beam are denoted by wq and wcr, re- 
spectively, where the plasma frequency of a plasma with 
species i is defined by 



E 



ATrqf 



dV 

lip') 



Mp') 



(1) 



where qi,mi, fi{p) are the species' charge, mass and mo- 
mentum distribution. It is assumed that ujq ^ cjcr, and 
that the CRs have a small but finite spread in the velocity 
directions. 

In this section we analyze the linear growth of unstable 
modes with long wavelengths, k -^ ujq. We start by 
considering a beam with no transverse velocity spread in 
§ 12.11 We show that the entire k-space regime considered 
is unstable and provide a simple analytic expression of 
the instability growth rate. The effects of a spread in the 
velocity directions of the CRs are discussed in 

2.1. No spread 



Consider the simplest case in which all the particles in 
the beam propagate in the same direction and are ultra 
relativistic (with an arbitrary energy distribution). It 
is straight forward to write the full dispersion equation 
which turns out to be a sixth order polynomial equation 
for oj with real coefficients (cf. § El and e.g. lAkhiezen 
119751 § 6.4.2). Four of the six solutions for w are small 
perturbations, of order Wqj^/wq, of the four cold plasma 
oscillating modes uj — ±wo, ±(wo + fc^)°'^, and are real 
(stable). The two remaining solutions have a non zero 
imaginary part and therefore are complex conjugates of 
each other. Hence for each k there is one unstable mode. 
Below we derive, directly from the Maxwell- Vlasov equa- 
tions, an approximate expression for the growth rate of 
this mode, Eqs. (fTO)) . and (|lll) . More details and a nu- 
merical solution for the dispersion equation are given in 

§E1 

For any axi-symmetric distribution of particles the lin- 
ear modes can be separated into modes having an electric 
field in the x — k plane, where x is the axis of symme- 
try (and magnetic field perpendicular to this plane) , and 
modes with an electric field perpendicular to x and k 
(cf. §E]). The unstable mode has an electric field, E, in 
the X — k plane, and a magnetic field, _B, perpendicular 
to this plane. The electrical currents carried by the cold 
plasma and the CRs, as derived by the Vlasov equations, 
are respectively given by 



47rJn = ^E, 



and 



47r JcR,± = «w^r(-B - Ei_)/n- 



AnJc 



R:|| 



^CR 



{B-E_L)k±/n^ 



(2) 



(3) 
(4) 



where subscripts _L and 1 1 correspond to components that 
are perpendicular to the beam and parallel to the beam 
respectively, fl = ku — uj, and where we used the ultra 
relativistic approximation /3 = 1 for the CRs. 

By neglecting the displacement current dtE, compared 
to the current carried by the cold plasma [using cj ^ cjq 
which is self consistently implied by the result, Eqs. pT " 
and (fTTj) ]. the Maxwell equations read 



ik^^B = 47r(JcR,|| +-^"0.11); 
-ik\\B = 47r(JcR,_L + Jq,±)'-, 
iojB = —ik±E\\ + ik\\E_i_. 

Equation ([7]) can be written as 



iuj{B~E±_) = -~ik^E\ 



II' 



(5) 
(6) 

(7) 



(8) 



by self consistently neglecting VtE± as follows. In regimes 
where f2 <C fcii, this term is negligible compared to the 
term k\\E^. Otherwise, where fJ > A:||, Eq. pTJ) . implies 
that U. ^ Lo and Eqs. (O and Q imply that E^ -^ B, 
making the term QE±^ negligible compared to uB. 

By neglecting the term ik^B in Eq. ^, compared to 
JcK,\\ [using Eqs. (|4|) and (fT0|l ], and substituting for 
•^CR,||i •^Clh^lh equations (H]), ([2]), and ^ respectively, 
equation ^ becomes 



WcR-Sfc-L/^^ = -uj'oB/kj_. 



(9) 




Fig. 1. — Growth rate of unstable modes for a delta function 
momentum distribution of beam particles with plasma frequency 
'^Cr/'^O = 0.01 and Lorentz factor 7 = 5000. The solid lines are 
the exact solutions of the dispersion relation, evaluated numerically 
(cf. § |A1) The dashed line (covering the solid line with fen = 
O.OSiJo) is the analytical approximation given in equation l|10|l . The 
curvature at low values of fc 1 is due to the electrostatic mode. 



The solution of this equation for u is 



Im{uj} = ± ■fcj_ = ±T]Q{k±); 



Re{uj} = k 



OJQ 



II- 



(10) 

(11) 

Eq. ([5]) determines the ratio between E± an d B. The 
resul t, Eqs. (fTO)) a nd (fTTI). agrees with that o f (jAkhiezeii 
[T975I § 6.4.2) and (|Lemoine fc PeUetieii [20091 ) in the rel- 
evant k space regimes. 

By retaining the term ikj^B in Eq. ([5]), and assum- 
ing finite Lorentz factors, the dispersion relation can be 
similarly solved resulting in 



Im{ui} = ±iLJcRi 



^^^^ 



where averaging is over qf fi{p) / {rriij) [cf. Eq. 
solution, which is valid for all k with ku <^ 1, 
Eq. pU|) in the regime 



(12) 



©]. This 
reduces to 



fell (7^ 



-) < fc_L < Wq. 



(13) 



In figure ([T]) the analytical approximation for the 
growth rate given in equation (|10|) is compared with a 
numerical solution of equation (IA2I) for different modes. 



For illustration modes with k\ 



II 



Wo that are not ana- 
lyzed in this paper are shown. As can be seen in the 
figure, the analytic expression provides an excellent ap- 
proximation in the relevant k-space regime. 

2.2. With spread 

For a velocity distribution of the beam particles that 
differs from a delta function, Eq. ^ should be replaced 
with 

ikj_B = -iuj't^YiB ki_ { — , ^ ,^ ) - iujlB/kj_, (14) 



{n - k^p^f 



or 



1 



(fi - k^p^_f 



1 



(15) 



where the same approximations leading to Eq. (jH]) were 
used, fl = fc||/?|| — w, and r^g is the growth rate for a delta 
function momentum distribution Eq. (jlOp . Averaging is 
carried over the velocity distribution function. 



1 



(f]-A;x/3±)2 
where 



4r i-00 (^ - k^p^y 



2^/3^, (16) 

/(/3±) = E ^ / ^,Mp'W'^ Pi-)- (17) 



It is to be understood that whenever a singularity is en- 
countered, the expression should be evaluated at J7 — > 
J7 — ie in the limit e — >■ 0"*". 

Consider first the following ID rectangular distribu- 
tion: /3, = /3|| = /3cos(_0), py = p^= /3sin(^), /3, = 0, 
w ith 6 un ifor mly distributed between ±A0, and j3 = 
^\ — 7~2 (see lAchterberg et al.l 120071 for a discussion 
of a distribution with two identical counter streaming 
beams). For small angular spread the growth rate can be 
approximated analytically by neglecting the variations of 
Px (assuming fix — P)- Under this assumption Eq. (fT6l) 
reads 

' ^ ' (18) 



(O - k^pi_f 



the dispersion relation [Eq. 
1 



(fci/3A6i)2 - VL 



2 ' 



{k^l3Mf 
and the solution for Vl^ is 



1^2 



reads 
1 



(19) 



n^ = {k^pMf 



vl- 



(20) 

For illustration, the growth rates of the unstable modes 
are shown in figure ([2]) as a function of the spread A6' for 
the ID rectangular distribution considered above. The 
solid lines in this figure are the results of a semi-analytical 
calculation in which the velocity integrals where evalu- 
ated analytically, and a continuous solution for the dis- 
persion equation as a function of the spread was found 
numerically, starting from the solution for a delta func- 
tion distribution. The obtained solution is verified to be 
the fastest growing one at a given k, by comparing it 
to the results of a '1-mode' PIC simulation (shown as x 
signs, cf. §[B|. These results are compared to the growth 
rate estimate of Eq. (|20| (dots) and to the growth rate 
of the delta function 770 distribution (dotted line). As 
can be seen in the figure the estimate given in Eq. (ppj) 
is in very good agreement with both the simulation and 
the semi-analytical calculation. 

Note the following features of the solution given by Eq. 
(123: 



For a sufficiently small spread, /S.6 < A^crit there 
exists an unstable mode with Re{u)} = fcii/?, where 
the condition for instability is 



1, 



k^i3Ae < 770, 

I3A9 < wcr/wo 



(21) 
(22) 



2. The growth rate, Im{Q) monotonically decreases 
from 770 to as the spread A6 is increased. 



x10" 



or 



-1/2 




k^/cOg= 0.4 



10" 



Fig. 2. — Growth rate calculations in the presence of a spread in 
the transverse velocities of the b eam particles. Different methods, 
as described in the text in § 12.21 were used to calculate the growth 
rates. Results are given for {fci|/a;o = 0.05,kj_/u)Q = 0.4(0.2)} 
Weibel modes in red (blue) of a beam with 7^ = 5000, ujcr/i^o = 
10~^ and different spreads, A9. Solid lines, dots and x symbols give 
the g row th rates obtained using the semi-analytical calculation, 
Eq. II2OI 1. and a '1 mode PIC simulation' (cf. § [BJ respectively. 
Dotted lines show the growth rate for a delta function momentum 
distribution (no spread, A0 = 0). 



3. For A9 > A^ciit the mode becomes stable. 

The criterion for instabihty given in (1211) can be inter- 
preted as the simple requirement (Akhiezer 1975i , § 6.4.2) 
that the particles in the beam do not move in the direc- 
tion perpendicular to the beam a distance exceeding the 
"wavelength of the mode during one e-folding time (of the 
delta- function instability) . 

We next show that most of the features of equation (|20p 
are generic to a large class of axi-symmetric distributions 
of CRs with a small spread. 

For purely imaginary O = —irj, Eq. (|16p reads 



1 



in-k±p±y 



kipi 



T 



iv' 



klPlf 



f(P±)dPi 



and has a zero imaginary part for all real 
limits 1] —i' 0+,oo the term {{ft — k±P±)~ 
—kj^ \P^) and respectively, where 






dPi 



PI 



dlii 



(23) 

In the 

roes to 



(24) 
(25) 



The last equality follows from the fact that 

7 r^dx = 0. (26) 

The following features of Eq. ([T5|) generalize the fea- 
tures of the ID distribution considered above with the 
condition for instability, 



generalizing the condition (PT|) . 



(28) 



1/2 



k^{Pl'y'''<m, 



(27) 



1. For a small spread, {P^) < rjo/k^ (equiv- 
alent to kj^ {P'l') ^ l/'^o)' there exists an un- 
stable mode, with Re{Lo} = kufH. To see this, 
note that the term ((J7 — fcj_/3_L)~^) continuously 
changes from — fc^ (/3~^) to as 77 changes from 
to 00, and must be equal to — ?7q for some positive 
■q. 

2. Assuming that the term {{ft — A:j_/3_l)^^) is mono- 
tonically .increasing with 77, for a marginal spread 
{(3J^) -^ Vo/k± the mode becomes stable 
?7 ^ 0. This suggests that for larger spreads the 
mode is stable. 

As for the 1-D distribution considered above, the cri- 
terion for instability given in (|27|) can roughly be inter- 
preted as the requirement that the particles in the beam 
do not move in the direction perpendicular to the beam 
a distance exceeding the wavelength of the mode dur- 
ing one e-folding time. Note however that the simple 
expression, A6f3, in the 1-D case, which is equal to the 
maximal velocity of the particles in the direction per- 
pendicular to the, beam, is replaced by the non trivial 
average, {P±'^) [the velocity average defined in ([M]) ]. 

which has a less obvious meaning. 

3. APPLICATION TO COLLISIONLESS SHOCKS 

We next discuss the possible application of the above 
results to the study of long wavelength magnetic field 
generation in the far upstream of collisionless shocks. 
The application of linear analysis of homogenous distri- 
butions to the non-homogenous and non-linear problem 
of collisionless shocks is far from being trivial. Further- 
more, the long wavelength modes studied above are not 
the fastest growing modes, and can be affected by the 
faster growing, short wavelength modes once the latter 
reach the non-linear stages. Nevertheless, the analysis of 
linear growth of long- wave length modes is an important 
step in the study of long- wave length magnetic field gen- 
eration and can be used as a basis for comparison once 
more accurate calculations are made (e.g. PIC simula- 
tions). In addition it is possible that some of the main 
features of the linear modes also appear in the more com- 
plicated shock scenario. 

Consider a shock with Lorentz factor Fsh ^ 1 propa- 
gating into a cold plasma with particle density rip and 
plasma frequency wq. We first assume that all parti- 
cles in the plasma have the same mass (e.g. electrons- 
positrons) and then generalize the results to an electron 
proton plasma. 

Assume that high energy, shock accelerated cosmic 
rays carry a fraction ~ tp of the post-shock energy. 
In the shock frame, the cosmic rays are not highly 
beamed and have an energy distribution ncr^si> 7s) ^ 
epTsh'n'o{'js/^sh)~^^^ with p ~ 2, where 7s is the Lorenz 
factor of the cosmic rays in the shock frame. 

We wish to analyze the generation of long wavelength 
magnetic fields in a region surrounding a point, x, deep 
in the upstream due to the interaction of the CRs that 
reach this point and the incoming upstream particles. We 



assume that this point is reached by a substantial fraction 
of the CRs that have Lorenz factors larger than a space 
dependent minimum 7s (a:), and study the instabilities in 
a simplified homogenous model of the upstream frame. 

The cosmic rays in the upstream frame are beamed to 
an angular separation of ^ l/Tsh, and have a plasma 
frequency 



^CR ~ 6l/2(^^/r,,,)-P/2a;o, 



(29) 



where ujq is the upstream plasma frequency and different 
values of 7s represent different positions in the simplified 
picturcQ. Deep in the upstream, where 7s (a^) ^ Tgh, 
the plasma frequency of the cosmic rays is much smaller 
than that of the upstream, and the analysis of section § [2] 
holds. 

Eqs. (fTO|) and ([29]) imply that long wavelength modes 
will grow with a growth rate of approximately 

m = k^el^\-fs/T,hrP/\ (30) 

It is useful to compare the e-folding time 77(^ 
of this instability to the time it takes the ambi- 
ent magnetic field Bq to deflect a cosmic ray par- 
ticle by an angle of l/Fgh back to the downstream, 
TRfi '^ Tshlsmc/ieBoTsh) = ^smc/{eBo). The ratio of 
the two times is 



Tr,,uVo ~ r. 



sh 



fc± 1/2 - 
Wo ' 



,^„/'(7s/^s.)l-p/^ (31) 



where cba ^ i3o/(87rnomc^) is the ratio of the energy 
density in the ambient magnetic field to the upstream 
rest mass energy density. For small values of the mag- 
netic field in the upstream, there will be a large range of 
k vectors 

k. > T-,%^''4l (32) 

for which the instability will grow on time scales that are 
much shorter than the deflection time of the particles. 

Far away from the shock in the upstream, where wcr is 
sufficiently small, the instability will be suppressed due 
to the \/T sh spread in the cosmic rays transverse veloc- 
ities. Using Eqs. (gT]) [or dH])] and dSO]), the modes are 
unstable only at locations in the upstream where 

ls{x) <Tsh(ey^Tsh)'\ (33) 

For illustration, the growth rate of a specific unsta- 
ble mode is shown in figure ([3]) as a function of 7s. For 
simplicity, the momentum distribution is assumed to be 
the same as in § 12.21 with parameters chosen in terms of 
shock parameters as: ±A0 = ±.1/Tsh (with Fgh = 25), 
upstream Lorentz factor 7„ — 2'^s^shi and beam plasma 
frequency, ujcr as given in Eq. ([2^ with p — 2, and 
0.1. In the figure we also show the results of the 



semi-analytical calculation (cf. § 12. 2p . the growth rate 
estimate of Eq. pO| (dots) and the results of the 1- 
mode PIC simulation (cf. § El x signs). For compari- 
son the growth rate tjq of the delta function distribution 
(cf. § 12.11 dotted line), and the maximal 7s from equa- 
tion ([55)1 (red circle, y axis value arbitrarily set to 0) 
are also shown. As can be seen in the figure the delta 
function momentum distribution result is a good approx- 
imation at small angular spread, and the angular spread 



^ Note that the energy carried by the cosmic rays in this frame 
of the upstream particle rest mass energy den- 
eg = riQmc^ . 



greatly exceeds that ( 
sity for Vs > f-J i £( 




'cr,s 

Fig. 3. — Growth rate for a specific Weibel mode, with A:||/a;o = 
0.05 and k_i_/uio = 0.5, shown as a function of the shock frame 
Lorentz factor, 7s, of the beam particles. The beam particles have a 
ID rectangular phase space distribution (cf. ij |2.2l l with parameters 
as discussed in the text in § [3] The solid line, dots and x symbols 
give th e gro wth ra tes o btained using the semi-analytical calculation 
(cf. §[2l2ll. Eq. ^, and a '1 mode PIC simulation' (cf. § [Bj 
respectively. The dotted line is the growth rate obtained for a delta 
function distribution, 770 . The red circle shows the value of 7CR,s 
given by equation II33I I. beyond which the modes are predicted to 
be stable. 



suppresses the growth rate at high cosmic ray Lorentz 
factors in accordance with the estimate in equation p3p . 
Next, consider a shock propagating into a electron- 
proton plasma. The equation for the plasma frequency 
of the beam ((29| should be replaced by 



WCR 



cl/2,- 



(me/mp)^/^(7s,p/rsh)-^/'wo 



P/2,, 



(34) 



where 7s,p is the minimal CR proton shock frame Lorentz 
factor at the position considered and where the contri- 
bution of the electron CRs at a given particle energy was 
neglected. Given the transverse velocity spread, \/Tsh, 
of the CRs Eqs. (HH]) and ^F^ imply that the long wave- 
length modes do not grow for moderate shock Lorentz 
factors (see also iLvu barskv & Eichlcr j2006i) 



Fsh < lOOe. 



1/2 
p, — 1' 



(35) 



where ep = 0.1ep__i. For higher shock Lorentz factors, 
the long wavelength modes grow at locations where the 
minimal shock frame CR energy, es{x) (equal for CR 
electrons and protons), is sufficiently low [cf. Eq. ((33)) ] 



^s(.x) < Tshrrip (ep/^Fsh 



2/p 



imJmp)'/P. (36) 



4. DISCUSSION 

In this paper the growth rates of the long wavelength 
unstable modes arising from the interaction of a beam of 
ultra-relativistic CRs and an unmagnetized cold plasma 
were calculated. We have shown that in the ultra- 
relativistic limit all long wavelength modes are unstable 
with a growth rate 770 = fc^wcR/wo [Eq. PH)) ]. as long as 
the spread in the transverse velocity distribution of the 
beam is sufficiently small. An extension of this result, 
for finite Lorentz factors and k± ^ wq, is given in Eq. 



6 



p2|) . For large transverse velocity spreads the instability 
is suppressed. The condition for instability was derived 
for a large class of velocity distributions [Eqs. (|28p and 

®]- . . . 

The possible application of these results to the interac- 
tion of CRs with the incoming plasma in the far upstream 
of unmagnetized collisionless shocks was addressed in §[31 
In a shock propagating into an electron positron plasma, 
the instability grows in upstream regions where the min- 
imal shock frame CR Lorentz factor is sufficiently low, 



Is 



< r- 



1/2 



ijCp [cf. Eq. ([33l) ]. Farther upstream of the 
shock the instability is suppressed due to the low density 
of the CRs. If these instabilities dominate the particle 
deflections responsible for particle acceleration, the accel- 
eration of CRs to shock frame Lorentz factors exceeding 
Tsh^p [cf. Eq. dSS])] may be obstructed. 



The long wavelength instabilities considered do not 
grow in the upstream of a shock propagating into an 
electron proton pla sma having Lo rentz factors < 100 (cf. 
Eq. ([35l) , see also iLvubarskv fc Eichlcr 2006) implying 
that: 1. The modes considered are not important 
in the relativistic shocks responsible for the observed 
GRBs afterglow emission (except possibly at the earliest 
stages, when Tgh ^ 100); 2. Particle acceleration in 
electron-positro n plasmas suggeste d by the results of 
PIC simulations 'Soitko vskvl (|2008bl ). where these modes 
can grow, may differ from acceleration in electron-proton 
plasmas. 
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APPENDIX 

A. FULL SOLUTION OF THE DISPERSION EQUATION 

Consider the case, described in §[2l where the velocity distribution of the particles is axi-symmetric (around x). In 
this case averages over the velocity distributions of the type ( \), and (/3^) , where _L, z are directions perpendicular to 
X, are zero. As a result, terms in the beam susceptibility (cf. iMelrosdl 19861 ). XCR containing such averages cancel out, 
and the susceptibility will only have non-diagonal term in the plane defined by beam direction and the wavenumber 
vector, k (the x — k plane). In this plane the beam susceptibility, the electro-magnetic field susceptibility, and the cold 
plasma susceptibility have matrix forms which are respectively 



2l3kn 



{k^ 



XCR 



^CR 



Pk±ui^ 

n 



TF" 




Xem 



j^-k] 
k\\k± 



k\\k±_ 

,)2 _ i.2 



Xo = Wq, 



(Al) 



where J7 is defined after equation (|3]), and the magnetic field B is perpendicular to the x — k axis. In the remaining 
axis only the plasma frequency term remains and this axis gives rise only to the plasma oscillation modes. 

In the absence of the beam the dispersion equation, det {xo + Xem} = 0, is a forth order polynomial equation in 
ui with real coefficients which has the following 4 solutions, u) = ±ojo,±{uJq + /c^)°-^, that represent respectively the 
plasma oscillation and the electro-magnetic mode. In the presence of the beam the dispersion equation becomes 



det {xo + XCR + xem} = 0, 



(A2) 



which is a 6'th order polynomial equation in uj with real coefficients. This equation is a slight perturbation of the 
original dispersion equation with 6 distinct solutions. Four of these solutions are slight deviation from the solutions 
of the original dispersion equation, while the other two solutions, for the k space regime that is considered above, are 
the solutions discussed in § 12.11 The full dispersion equation can be solved numerically for any k and these solutions 
are shown in figure ([l}. 

B. 1-MODE PIC SIMULATIONS 

The solution described above was numerically verified to be the fastest growing one at several wavevectors, k. This 
was done by performing an efficient PIC simulation in which only one k mode is treated while the rest of the modes are 
neglected. In this simulation, the electric field, magnetic field, and electric current are sinusoidal with a given k value 
and a time dependent amplitude, while the cosmic rays are treated as particles with continuous position (ID along k) 
and momenta (3D). The upstream, which is assumed to have a delta function momentum distribution, is written in 
terms of fluid quantities in the linear approximation and is likewise sinusoidal with the given k. The time dependent 
amplitude of the electric current carried by the cosmic rays is derived from the distribution of the particles by 



where gj,/3., and x^ are the charge velocity and position of the particle j. As such, the simulation is only accurate for 
the linear regime (the amplitude at which the mode becomes non-linear can in principle be identified) . The growth rate 
is calculated by fitting the time evolution of the amplitude of the mode with an exponential function in time. With 
this method only the fastest growing mode at the given k is accounted for. This is a direct, physically transparent 
method of studying the linear regime of a single wavevector, with no restriction on the 3D velocity distribution, which 
is time and memory efficient. 

For illustration, the results of the simulation used for figure ([3]) are presented in figure (jBip . As an underlaying 
quantity for the fitting we use the amplitude of the electric current, jk, which is normalized to the current jo = gngc. 




Fig. B1. — The growth of a {kii/uiQ = 0.05, A:x/'*^o = 0-5} Weible mode in a 1-mode PIC simulation described in the text in §(5] The 
parameters chosen for this simulation were taken to reproduce the 7s = 25 (ijJq-h = 0.1, 711 = 1250) point given in figure ||3|. The green 
curve is the logarithm of the amplitude of the current vector, j'fc, obtained from the simulation, and the blue line is a linear fit to this curve 
in the range 19.6 < tuio < 44.5. Its slope, 0.131, represents the growth rate. 



where ng is the upstream density, and q is the electron charge. As can be seen in the figure, the current jk grows 
exponentially with time and the growth rate of the instability is easily obtained from the fit. The accuracy of the growth 
rate obtained numerically is high despite the fact that only 10^ particles (cosmic rays) were used. The saturation level 
of jk can also be obtained from this figure. However, since all quantities, except those related to the cosmic rays, are 
treated in the linear regime, this saturation is representative only if in reality the saturation is governed by the cosmic 
rays. 
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